A class of well-posed approximations for constrained second order hyperbolic equations
نویسنده
چکیده
The purpose of this paper is to present a new family of numerical methods for the approximation of second order hyperbolic partial differential equations submitted to a convex constraint on the solution. The principle is a singular modification of the mass matrix obtained by the mean of different discretizations of the solution and of its time derivative. The major interest of these methods is that the semi-discretized problem is well-posed and energy conserving. Numerical experiments show that this is a crucial property to build stable numerical schemes.
منابع مشابه
ar X iv : g r - qc / 0 10 60 85 v 1 2 6 Ju n 20 01 Some Mathematical And Numerical Questions Connected With First And Second Order Time Dependent Systems Of Partial Differential Equations
There is a tendency to write the equations of general relativity as a first order symmetric system of time dependent partial differential equations. However, for numerical reasons, it might be advantageous to use a second order formulation like one obtained from the ADM equations. Unfortunately, the type of the ADM equations is not well understood and therefore we shall discuss, in the next sec...
متن کاملar X iv : g r - qc / 0 10 60 85 v 2 7 J ul 2 00 1 Some Mathematical And Numerical Questions Connected With First And Second Order Time Dependent Systems Of Partial Differential Equations
There is a tendency to write the equations of general relativity as a first order symmetric system of time dependent partial differential equations. However, for numerical reasons, it might be advantageous to use a second order formulation like one obtained from the ADM equations. Unfortunately, the type of the ADM equations is not well understood and therefore we shall discuss, in the next sec...
متن کاملNumerical stability for finite difference approximations of Einstein's equations
We extend the notion of numerical stability of finite difference approximations to include hyperbolic systems that are first order in time and second order in space, such as those that appear in Numerical Relativity. By analyzing the symbol of the second order system, we obtain necessary and sufficient conditions for stability in a discrete norm containing one-sided difference operators. We pro...
متن کاملError analysis of implicit Euler methods for quasilinear hyperbolic evolution equations
In this paper we study the convergence of the semi-implicit and the implicit Euler methods for the time integration of abstract, quasilinear hyperbolic evolution equations. The analytical framework considered here includes certain quasilinear Maxwell’s and wave equations as special cases. Our analysis shows that the Euler approximations are well-posed and convergent of order one. The techniques...
متن کاملPerfectly Matched Layers for Second Order Wave Equations
Numerical simulation of propagating waves in unbounded spatial domains is a challenge common to many branches of engineering and applied mathematics. Perfectly matched layers (PML) are a novel technique for simulating the absorption of waves in open domains. The equations modeling the dynamics of phenomena of interest are usually posed as differential equations (or integral equations) which mus...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 2008